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We demonstrate sum-frequency generation in a nonlinear whispering gallery mode resonator be- 
tween a telecom wavelength and the Rb D2 line, achieved through natural phase matching. Due to 
the strong optical field confinement and ultra high Q of the cavity, we achieve a 1000-fold enhance- 
ment in the conversion efficiency compared to existing waveguide-based devices. The experimental 
data are in agreement with the nonlinear dynamics and phase matching theory in the spherical 
geometry employed. The experimental and theoretical results point to a new platform to manipu- 
late the color and quantum states of light waves toward applications such as atomic memory based 
quantum networking and logic operations with optical signals. 



Unique properties of optical whispering gallery mode 
(WGM) resonators have been attracting a continuously 
increasing attention in the recent years. An extremely 
high Q-factor and narrow linewidth of such resonators, 
in combination with their excellent mechanical stability, 
gave rise to a number of optical filters and laser stabi- 
lization applications [3-|g|. On the other hand, the solid 
composition of WGM resonators provides sensitivity to 
the environmental changes that affect the solid's opti- 
cal properties, or perturb the evanescent optical field 
surrounding the resonator, and enables various sensor 
applications [2H3I- Among other perturbations, WGM 
mechanical oscillations strongly coupled with the opti- 
cal fields may open up exciting possibilities in cavity- 
optomechanics with strong photon-phonon coupling. 

For WGM resonators possessing nonlinear optical sus- 
ceptibilities, very efficient nonlinear optical processes are 
enabled at pump powers orders of magnitude smaller 
than are required in traveling- wave devices. Thus far, 
a number of second- and third-order nonlinear processes 
have been demonstrated in these devices including four- 
wave mixing |10| - KL4| , frequency doubling [15l - [l7| , tripling 
[18l IHI. quadrupling |2(| an d parametric down conver- 
sion [Hriiji . 

In this Letter, we report sum-frequency generation 
(SFG), where one coherently up-converts the optical sig- 
nal frequency with a local oscillator, in WGM resonators. 
SFG can be employed for efficient room-temperature 
detection of far- infrared, and even sub-THz, photons 
[2a . l2q . Such conversion has been shown to preserve the 
quantum state of the input signal [2|| [27|. It also can 
lead to efficient manipulation of the color and shape of 
single-photon signals |28j |. Furthermore, the discrete and 
sharp nature of the resonance lines in such devices can 
greatly suppress incoupling of Raman noise, and poten- 
tially lead to new optical tools for mode discrimination 
and reshaping of narrowband quantum signals [29| . We 
demonstrate, for the first time, highly-efficient, naturally- 
phase matched SFG between a telecom wavelength (1560 



nm) and the Rb D2 line (780 nm) at pump power lev- 
els in the 100-//W range, considerably smaller compared 
to waveguides due to the enhanced interaction length in 
these devices [30j. Low noise and operation power, along 
with compactness and stability, make these devices suit- 
able for applications such as logic operations j3lM33| and 
networked quantum communications using atomic mem- 
ories [3J,|35(. Given these benefits, we decided to verify 
the WGM SFG in an experiment illustrated in Fig. [TJ 
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FIG. 1: Experimental setup. PD: Photodetector, M: Mir- 
ror; DM: Dichroic Mirror; PBS: Polarization Beam Splitter; 
HWP: Half- Wave Plate; DAQ: Data Acquisition Unit. 10X 
Objective Lenses were used to focus the lasers onto the prism- 
resonator interface and to collect the output light. 

We observed sum-frequency generation in a MgO- 
doped Lithium Niobate z-cut microdisk (R ~ 0.6mm) 
evanescently coupled to a diamond prism. Diamond- 
polishing was used to obtain Q > 2 x 10 7 and Q > 4 x 10 7 
for the signal and pump waves, respectively. A 780nm 
narrow (below 300 kHz) linewidth tunable-diode laser 
was used as the signal and a DFB laser provided 1560nm 
pump. Two input waves are ordinarily polarized to 
achieve the Type-I phase-matching whereas the upcon- 
vertcd wave at 520nm is detected in the extraordinarily 



polarization. The two input waves were combined on a 
dichroic mirror and focused onto the prism-resonator in- 
terface by an objective lens. A lateral offset of the beams 
before the lens allowed for optimizing of the pump and 
probe in-coupling angles individually. A similar lens was 
used to collect the output light, with a dichroic mirror 
separating the pump and the signal while a polariza- 
tion beam splitter separating the sum-frequency wave. 
All three optical powers were measured by photodctcc- 
tors, whose signals were fed into a data acquisition unit. 
As the laser frequencies were continually swept at 50Hz 
across several linewidths, the signal and pump WGMs 
were tracked by a software program continuously adjust- 
ing the lasers central wavelength to follow their respective 
WGMs. The program ensured that the pump and signal 
WGMs were in the center of the sweeps and that they 
were pumped simultaneously In addition, the top of the 
resonator was coated with silver paste and temperature- 
controlled to allow for electro-optic and thermal tuning 
for the SFG phase matching. 

Having the phase matching achieved, we measured the 
sum-frequency output power for various input pump and 
signal powers. In this measurement we used a slowly 
varying electro-optic bias voltage to record the transient 
peak SFG efficiencies circumventing the relative instabil- 
ity of the temperature controller. For the same reason 
each data point represents the average of three consecu- 
tive measurements. 

The signal and pump waves were critically coupled and 
over-coupled, respectively. The longer wavelength pump 
wave was coupled stronger than the signal due to the na- 
ture of the evanescent coupling. Also, due to the spatial- 
mode mismatch between the input Gaussian beam and 
the WGM profile, in this measurement we achieved a 
critically-coupled contrast of 48%. We took this into ac- 
count by only utilizing the in-coupled powers for our the- 
oretical analysis. 

We observed efficient sum-frequency generation with a 
maximum in-coupled pump power of only 1.22 mW. In 
Fig.[2]we plot the out-coupled SFG efficiency As we var- 
ied the signal powers, saturation of the peak conversion 
efficiency was observed in all cases with sub-milliwatt 
incoupled powers, instead of a cyclic behavior that is 
observed in the traveling- wave configuration |36[. Sim- 
ilar saturation has been observed in frequency-doubling 
WGM experiments [la. 1 161] . At higher pump powers, an 
additional nonlinear loss for the signal wave is created 
due to its upconversion, leading to the reduction of its 
internal Q-factor and the coupling contrast. As a result, 
a smaller portion of the signal wave enters the resonator 
and the SFG efficiency is reduced. This behavior is a 
manifestation of the "coherent" quantum Zeno effect for 
the signal wave, where the "potential" for the upconver- 
sion decouples the signal field from the cavity [321 ] . 

A theoretical description of the SFG in WGMRs is 
warranted because manifestly different behavior is ob- 
served at high pump powers compared to a traditional 
traveling-wave geometry Neglecting Rayleigh backscat- 
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FIG. 2: The out-coupled SF emission is measured and nor- 
malized to the various input signal powers. Symbols represent 
the experimentally measured data, and the solid lines are the- 
oretical fits. Note that the theory predictions for 2.5{iW and 
5fiW are identical because the nonlinear loss is negligible at 
such powers. 



tcring and assuming linearly polarized fields, we find the 
scalar equations of motion in the cavity are: 
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where Q c and Q l denote the coupling and intrinsic Q- 
= (5o^ + 5o^ -z A") for M = s,p,/ indicating 



factors, 



respectively signal, pump and the sum-frequency, and 
Q = 7^31 / dVEfE*E* is the internal conversion effi- 
ciency of the SFG process. The input field operators, a M , 
are related to the output fields by 6 M = VTa^ + i.l-MrC^. 

Using quasi-static analysis we solve Eqns. (l)-(3) for the 
out-coupled sum-frequency field, |6/| 2 = \i./^-^-c s c p \ 2 . 

Before we provide the solution, a few salient points are 
to be noted. In this formulation, tt is the internal con- 
version efficiency of the SFG process, whereas only the 
SF out- coupled power is measured experimentally Since 
Q c is determined by the distance from the prism to the 
resonator d, we need only to fit to f2 and the intrisic 
Q- factor for the SF [37|. Moreover, these equations are 
difficult to solve in general for c s , c p , and c/ analytically. 
To guide us, however, we make the undepleted pump ap- 
proximation and then use numerical methods to acquire 
the generic solution. This approximation yields: 
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Note that the expressions for c p and c s are asymmetric 
due to the nature of the undepleted pump approxima- 
tion. These solutions indicate that we do not observe an 
oscillatory behavior in the frequency-conversion dynam- 
ics. For low pump and signal energies, the SFG output 



behaves linearly, \bfY 



ex \a v 



(see Fig. [2]), whereas 



at higher energies, the dynamics arc different, i.e., the 
upconversion and subsequent down conversion processes 
are asymmetric in this geometry [321 ] - 

Using the in-coupled pump and signal powers, mea- 
sured Q M and calculated Q\,Q c p for \i = s,p, we can 
provide a theoretical fit for the measured data with the 
fitting parameters £7 and Q\. Fitting our efficiency mea- 
surements at 2.5/iW signal inputs (where the undepleted 
pump approximation is valid) with Eq. ([5]), we can esti- 
mate J7, which does not vary as the same WGM triplet 
is studied throughout the experiment. Throughout this 
procedure, wc assumed that the SFG peak occurred 
when the three waves were exactly on resonance, i.e., 
A s = A p = A/ = 0. 

Wc then solved for Eqs. (l)-(3) numerically, which pro- 
vided the efficiency curves for data where the undepleted 
pump approximation is invalid. The total Q-factors for 
the pump and signal waves were not constant through 
the experiment, which wc account for by using Q\ as 
a fitting parameter in Fig. [2J Wc attribute this varia- 
tion in part to the unaccounted photorefractive effects in 
the resonator, caused by the SF. These effects were most 
evident for the highest signal power used, P s = 60[iW. 
Indeed, this particular data set in Fig.[5]also presents the 
largest data scatter and the worst agreement with theory. 

In spite of these theoretically unaccounted background 
processes, by using just the O-parameter, which deter- 
mines the shape of the efficiency curves, we are able to 
accurately model the signal and pump transmission spec- 
tra as well as the SF emission spectrum [33]. Moreover, 
we were able to calculate [37[ the fundamental channels 
modes overlap to give fltheor = 253 kHz, whereas the 
empirical value gives us £l exp t = 5 kHz. In calculating 
this latter value, we used Q\ = 3.25 x 10 7 , the intrinsic 
sum-frequency Q-factor and that it was strongly under- 
coupled to the resonator, -=rf = 164, see [371 ]- 

Not every pair of signal and pump WGMs can generate 
sum-frequency. SFG in a triple-resonant system requires 
the phase matching between these modes, which can be 
viewed as conservation of the integrals of motion deter- 
mined by the system's symmetry. 

Usually WGMs have nearly perfect spherical symme- 
try, so their eigenfunctions inside the resonator can be 



where r,6,ip arc spherical coordinates, L,m,q are az- 
imuthal, polar and radial mode numbers, respectively, 
Y™(9, cp) is a spherical harmonic and JL(k q r) is a spheri- 
cal Bcsscl function. The radial wave number k q is deter- 
mined from the boundary conditions. 

In [371 we show that many combinations of WGM 
triplets (called channels in the following) may lead to 
SFG albeit with different conversion efficiencies, whereas 
most efficient channel is the one that couples the fun- 
damental modes, i.e. such that q p = q s = qj = 1, 
L p — m p = L s — m s = Lf — irif = 0. 

Achieving a triple resonance for a selected channel re- 
quires tuning of each WGM's frequency such that the 
energy conservation ui p {L p ,m p ,q p ) + u> s {L Sl m s ,q s ) = 
ujf(Lf,m s + m p ,qf) is fulfilled to better than a WGM 
lincwidth. In lithium niobate resonators this can be 
achieved for ordinary pump and signal and extraordinary 
sum-frequency WGMs due to different temperature de- 
pendencies of the ordinary and extraordinary refraction 
indices. We find the phase matching temperatures [371 ] by 
iteratively solving the energy conservation condition us- 
ing the WGM dispersion equation 1391] and temperature- 
dependent Scllmcier equations [40| . Let us point out 
that for sub-mm resonators the resulting temperatures 
may significantly (by tens of degrees) differ from the bulk 
phase matching temperature, due to the geometrical, or 
waveguide, part of the WGM dispersion. 

In the experiment, identifying a WGM's q may present 
a considerable difficulty. Fortunately, the WGM's free 
spectral range (FSR) depends on its q much stronger than 
on L and m, because q affects the effective length of the 
resonator. Therefore a WGM's q can be inferred from 
the FSR measurements. We carried out the FSR mea- 
surements for the pump laser by frequency-modulation 
technique [41|. The results of such measurements car- 
ried out with 11 best-coupled modes within one FSR are 
shown in Fig. [3J The theoretical FSR values shown in 
Fig . [3] were derived from the WGM dispersion equation 
[391 ] . We fit the theory value for q = 1 to the smallest 
measured FSR of 32.362 GHz by varying the resonator 
radius from the initially measured 0.65 ± 0.01 mm to 
0.6505 mm. No fitting was done to match other q values. 

The theoretical FSR value for q = 1 WGMs at the sig- 
nal wavelength is 31.049 GHz. We found a high-contrast 
mode with a very close FSR value of 31.04 ± 0.004 GHz 
[43 ] . Coupling the pump and signal lasers to these WGMs 
we slowly varied the resonator temperature while moni- 
toring the SFG signal, and acquired the data for Fig. [2] 

We also found the neighboring SFG channels by tun- 
ing the pump and signal lasers across an integer number 
of FSRs. The wavelengths and phase matching temper- 
atures for these channels are shown in Fig. [4] together 
with the numeric simulation for the fundamental chan- 
nel. A good agreement was achieved by using the MgO 
concentration as a fitting parameter in the simulation. 
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FIG. 3: Theoretical FSR values for q = 1...8 WGMs at 1560 
nm (vertical lines), and the measurement results. 



This parameter affects the phase matching temperature 
by entering Lithium Niobate dispersion |40[. Unfortu- 
nately its value is not precisely known for our congruent 
wafer, except that it should slightly exceed the threshold 
value of approximately 5%. The fitting yielded a very 
plausible value of 5.63%. It should be pointed out that 
no other efficient SFG channels can fit the observations 
with any reasonable MgO concentration. However, many 
less efficient equatorial channels exist for the same tem- 
peratures [33]. The lower observed f2 suggests that the 
WGM triplet may not have been fundamental. 
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FIG. 4: Numeric simulation result for the fundamental 
SFG channel (black dots) and experimental observations (red 
stars). Projections emphasize a good agreement between the 
theory and experiment. 



To summarize, we have demonstrated triple-resonant 
SFG in a WGM resonator. We have extended the the- 
oretical analysis for finding the phase-matched WGMs 
inside the resonator and understood the nonlinear dy- 
namics of frequency conversion in the strong pump-signal 
coupling regime. The efficiency of this process in the res- 
onator is much higher than in a traveling- wave geometry, 
requiring sub-milliwatt powers for saturation. We expect 
to find applications of this interaction in the fields of spec- 
troscopy, optical communications and data processing, 
both at classical and at quantum levels. Other conceiv- 
able applications are for fundamental tests of quantum 
theory, e.g. in cavity optomcchanics and quantum non- 
demolition measurements. 
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I. ACQUIRING fi 



For low signal powers, the sum-frequency efficiency can 
be understood in the undepleted pump regimes (Fig. [TJ. 



are given by [lj for a TM WGM in a spherical resonator: 
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FIG. 1: The dotted line is the exact numerical solution for 
2.5/j.W signal in-coupled power, and the solid line is the unde- 
pleted pump approximation. Their good agreement enables 
us to estimate f2 and utilize it for higher signal powers, where 
the approximation is invalid and exact analytical quasi-static 
expressions do not exist. 



Since Eq. (5) in the main text is valid only in the un- 
depleted pump approximation, we can infer Q only for 
datasets where the signal powers are very small. Addi- 
tionally, we assume that the inherent nonlinearity is so 
weak that for both pump and signal powers P < 5(iW , 
there is no noticeable conversion (which is justified by all 
of our data sets, wherein the efficiency is nearly zero in 
this region) . Since we use the same WGM triplet for the 
entire experiment, and Q is an overlap measure between 
the excited WGMs, we operate under the condition that 
fl is the same for arbitrary input powers. However, what 
we can acquire from the measured data is Q and K/, 
where Kf = Wf(okr + 2Q°) - ^he coupling-Q factors, Q°, 



\x = s,p, where n s , n c are the indices of refraction in the 
resonator and prism respectively, R is the resonator ra- 
dius, 7 = y/k 2 (rig — 1), and d is the distance between the 
prism and resonator. We will use ([1]) as an order of mag- 
nitude estimate for our spheroidal resonator, and multi- 
ply it by nj, when computing Q c f for the extraordinarily 
polarized (TE) 520 nm sum-frequency wave, which is due 
to the different boundary conditions. Since we operate 
under critical coupling conditions, i.e., Q c s = Q\, for the 
X s = 780 nm signal wave, we can acquire an empirical 
value for d. In our experiment, for Q c s = Q\ = 4.48 x 10' 
leads to d = 70.5 nm. For this distance, for the 520 nm 
wave, the calculated QS = 5.34 x 10 9 . From data fitting 
we find Q\ = 3.25 x 10 7 , which leads to the undercou- 
pling factor Q C JQ\ = 164. This estimated intrinsic Q is 
very close to an earlier reported value of Q l — 4 x 10 7 
in the same material at a close wavelength. 

Therefore, the internally generated green light was cou- 
pled out very inefficiently. This leads us to believe that 
the photorefractive damage induced by the green light 
could have been significant. Moreover, from Eq. (1), 
we see that any changes in the index of refraction due 
to thermo-optic effects, light-induced charge transport, 
or even the electro-optic effect due to charge movement 
on the periphery of the resonator can significantly affect 
Q c „. The exponential dependence on d places stringent 
stability criteria, which may have been breached due to 
the mechanical hysteresis in the brass-oven that was em- 
ployed to allow for temperature-tuning, leading to a time- 
varying Q c . In fact, by studying the transmission curves 
over the duration of the experiment, we have seen that 
the total Q of the resonator does change. Therefore, we 
conclude that the f2 we acquire by assuming that Ay = 
and by also neglecting these concurrent effects, might be 
different slightly from the value acquired when these ef- 
fects are accounted for. This may further explain why 
the empirically measured internal conversion efficiency is 
lower by a factor of 50. 



II. OBSERVING AND FITTING THE WGM 
SPECTRA 

As discussed in the main part, our measurement pro- 
duced the signal and pump transmission spectra, and the 
SF emission spectra, such as shown in Fig. [21 Having esti- 
mated the tt, Kf parameters from the experimental data, 
we were able to acquire parameter-free fits for the trans- 
mission spectra of the resonator, also shown in Fig. [2] 
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FIG. 2: Dotted lines are the observed signal, pump and SF 
spectra, and solid lines are theoretical fits. 



III. WGMS OVERLAP AND SFG CHANNELS 
EFFICIENCY 



Spherical symmetry may be a good approximation for 
determining the eigenfunction shape (but not always the 
cigenfrequencies) in a spheroid WGM resonator whose 
aspect ratio is not too large, because the variation of the 
boundary from a true sphere in the region of significant 
WGM field is minute. Due to the spherical symmetry 
the phase matching in WGM resonators corresponds to 
conservation of photons' orbital momenta 0t5|- If £ 3> 1 
the overlap integral in D, breaks up into a product of 
the angular and radial part, former depending on L — 
m and latter on the q for all three modes. Each part 
also weakly depends on wavelengths and the three values 
of L. The angular parts give rise to the SFG selection 
rules corresponding to the Clebsch-Gordan coefficients. 
In particular, they enforce 



where N is an integer. The radial parts lead to no partic- 
ular selection rules, however as we will see they strongly 
favor the cases when q p + q s ps qf . 

For mm-size resonators the orbital numbers L,m are 
large, and evaluating the overlap integrals with Legendre 
polynomials and Bcsscl functions of such orders is im- 
practical. Appropriate asymptotic approximations need 
to be made. Let us introduce 



$(r,0,cp) PS ^ang{0)^rad(r)e l 



(5) 



inside the resonator, and ty(r,9,ip) = outside, thereby 
neglecting the evanescent field. The angular asymptotic 
for ([5]) is given by @ 

*ang(e) = -^=H L _ m (VL C Os(e)) exp (-1lcos(0) : 

V Z7T I ^ 



and radial asymptotic is 



(6) 



#, 



■<" - > 



L-k r (L + 1/2, q)r/R 



L -k r (L + 1/2, q) 

In (jfj|) Hh-m is the Hcrmitc polynomial of the order L 
m, and N a is the normalization factor, 



(<) 



n: 



H\_ m {yLx) exp {—La; 2 } dx. 



(8) 



In ||7J), R is the resonator radius, Ai is the Airy function, 
a q is its qth root (positive value), dimensionlcss wave 
number can be approximated [7[ as 

k r (L,q)^L + a q (L/2) 1 / 3 , (9) 

and the normalization factor is 

L -k r (L + 1/2, q)x 



n: 



R" 



Ai 2 



L - k r (L + 1/2, q) 



xdx. (10) 



The factors ((8|) and (|10|) provide the asymptotic wave- 
function ([3]) normalization: J dV\^(r, 9, ip)\ 2 = 1. 

Approximation ([5]) allows us to evaluate the radial and 
angular overlap integrals separately. These overlap fac- 
tors are shown in Fig. [3] in the normalized form. To 
determine the overlap integral for a particular channel 
one needs to take the appropriate angular part from 
the top of Fig. [3] and multiply it by the appropriate 
radial part from the bottom of Fig. [3] For the fun- 
damental modes, i.e. such that q p = q s = qj = 1, 
L p — m p = L s — m s = Lf — m/ = (the most ef- 
ficient SFG channel for our wavelengths and resonator 
size) \JdV^ f ^^* s \ 2 ps 2.66 x 10 6 cm~ 3 . For compar- 
ison, the fundamental WGM mode volume in our res- 
onator ranges from approximately 1 x 10~ 7 cm 3 for the 
SF to 4 x 10 -7 cm 3 for the pump. 



IV. PHASE MATCHING TEMPERATURE 



L/r> "T La 



m s = nif, 


(2) 


L s > L f , 


(3) 


L f = 2N, 


(4) 



Based on the previous section analysis, we can rank 
various SFG channels by their efficiency. We then nu- 
merically find the dependence of the frequency detuning 
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FIG. 3: The absolute-square of the angular (top) and radial 
(bottom) parts of the WGMs overlap integral are represented 
by the dot size as a function of the respective mode numbers. 
In both plots, the largest dot is normalized to unity. 



Au = uif — uj s — Lu p on temperature and determine the 
phase matching (Auj = 0) temperature as well as the 



WGM frequencies uif, ui s , uj p at that temperature, for 
each channel from the ranked list. 

In Fig. 3] we show the theoretical prediction for phase 
matching temperatures in the range between 120 and 150 
°C. The pump modes in this simulation have been as- 
sumed equatorial (L p — m p = L s — m a = 0) and with 
q s . p < 8. Both assumptions are justified by the high 
contrast of the observed WGMs and by the FSR mea- 
surements discussed in the main part of the paper. No 
assumptions have been made about possible SFG WGMs. 
From Fig. 2] we see that the most efficient [1, 1, 1] SFG 
channels have the phase matching temperature near 129 
°C, far from other channels of significant efficiency. How- 
ever, there are also many "minor" SFG channels at this 
temperature. In our experiment, we indeed observed sev- 
eral other SFG channels with the efficiency at least an 
order of magnitude below than the one we have been 
using. 



[2,4,5] -j— T - 
[2,3,4] J 
[1,3,3] 



%* * 



[1,2,2] 



[1,1,1] 



» «*» ao» 



» &&& }§ 9 '■• 




150 -TJ 

IX 
03 
(/> 
-145 CD 



.140 B" 

3; 

: 135 °=> 

i— t- 

CD 

3 

,130 "D 

r cd 
3 

125 C 

— \ 
CD 



780.4 

780.3 



-^^=^120 O 
1560.6"^ 



1559.7 



1560.0 



\Nave 



wnrfh ^ ] 



FIG. 4: Main equatorial SFG channels [q s ,q P ,qf] form dis- 
tinct "layers" in the temperature-wavelength space. Within 
each layer, the channels form the rows of constant L a and L v . 
Dot sizes represent the channels efficiency. 
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